Liquid crystal cells with "dirty" substrates 
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We explore liquid crystal order in a cell with a "dirty" substrate imposing a random surface pin- 
ning. Modeling such systems by a random-field xy model with surface heterogeneity, we find that 
orientational order in the three-dimensional system is marginally unstable to such surface pinning. 
We compute the Larkin length scale, and the corresponding surface and bulk distortions. On longer 
scales we calculate correlation functions using the functional renormalization group and matching 
methods, finding a universal logarithmic and double-logarithmic roughness in two and three dimen- 
sions, respectively. For a finite thickness cell, we explore the interplay of homogeneous-heterogeneous 
substrate pair and detail crossovers as a function of disorder strength and cell thickness. 
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Over the past several decades there has been consid- 
erable progress in understanding the phenomenology of 
ordered states subject to random heterogeneities, gener- 
ically present in real materials [l[ . While most of the fo- 
cus has justifiably been on the bulk heterogeneity, there 
are many realizations in which random pinning is con- 
fined to a surface of the sample. A technologically rel- 
evant example, illustrated as the inset of FigflJ is that 
of a liquid-crystal cell (e.g., of a laptop display), where 
a "dirty" substrate imposes random pinning, that com- 
petes with liquid-crystal ordering. Manifestations of het- 
erogeneous surface anchoring include common schlieren 
textures 0] and multistability effects Q observed in nc- 
matic cells, as well as long-scale smectic layer distortions 
in smectic cells 01 ■ 

In this Letter we study such surface randomly-pinned 
systems modeling them via a d-dimensional xy model 
with the heterogeneity confined to a [d — l)-dimensional 
surface. One might a priori expect surface pinning ef- 
fects to be vanishingly weak compared to the ordering 
tendency of the homogeneous bulk, and thus the xy order 
to be stable to weak surface disorder in any dimension. 

Our finding contrasts sharply with this intuition. 
Namely, our key qualitative observation is that the xy 
order on the surface of such (i-dimensional system with 
d < die = 3 is always destabilized by arbitrary weak 
random surface pinning 

Our finding can be understood from a generalization 
of the bulk Imry-Ma argument 0, 0] to the surface pin- 
ning problem. For an ordered region of size L, the in- 
teraction with the surface random field can lower the 
energy by E pin ~ V p VN p ~ A 1 / 2 L^/ 2 , where V p is 
a typical pinning strength with zero mean and variance 
A/ ~ ^p 2 /Co _1 (£o is the pinning correlation length) 
and N p is the number of surface pinning sites. Since 
surface distortions on scale L extend a distance L into 
the bulk, the corresponding elastic energy cost scales as 
E e ~ KL d ~ 2 , where K is the elastic stiffness. By com- 
paring these energies it is clear that for d < 3, on suffi- 
ciently long scales, L > £ L ~ {K 2 / A/)^ 3- ^ the sur- 



face heterogeneity dominates over the elastic energy, and 
thus on these long scales always destroys long-range xy 
order for arbitrarily weak surface pinning^]. A more de- 
tailed analysis extends the argument to 3D. The surface 
Larkin length scale beyond which the orientational or- 
der on the T = random surface (z = 0) is destroyed, is 
given by 

£ L ~ae cK2/A ?, ford = 3, (1) 

where c = 87r 3 is a nonuniversal constant, and a is a 
microscopic cutoff of order of a few nanometers in the 
context of liquid crystals, set by the molecular size. 

Although these (heterogeneous) substrate-induced dis- 
tortions from short scales x < £l only penetrate a 
distance £l into the bulk, decaying exponentially be- 
yond this scale, the bulk ordered state never recov- 
ers. That is, strikingly, we find that logarithmically 
rough distortions, persist an arbitrary long distance z 
from the dirty substrate, as summarized by C(x, z, z') — 
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FIG. 1: Average orientational order parameter ip(w, 0) (con- 
trolling light transmission through a liquid-crystal cell) at the 
random pinning substrate of a 3D Dirichlet cell of thickness 
w, illustrated in the inset. 
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(0(x, z) - 0(0, z')) 2 ), with (for x > £l in 2D) 



C 2Z /(x, z,z) ~ b 2 e 
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(62 is a weak function of 2z/£l), and in 3D given by a 
double- logarithmic asymptotics (fT2| . 

With an eye to liquid- crystal cell applications [3, 0|i we 
extend our analysis to a nematic cell of finite thickness 
w, with a heterogeneous front substrate, and a back sub- 
strate with a homogeneous Dirichlet (T>) or Neumann 
(A/ - ) boundary condition, that can be imposed by various 
substrate treatments, for example, surface polymer coat- 
ing and rubbing. As expected, for the D-cell long-range 
orientational order is stable to weak pinning, but with a 
value of surface orientational order parameter ip that ex- 
hibits a strong crossover as a function of w/£l as shown 
in Fig. (p}, that in 3D is given by 
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thin cell, a <C w <C £l, 
thick cell, w ^> £l, 



(3) 



where 7730 = 7r 2 /18 is a universal exponent and a = 
2tt 2 , rjL = 27r 2 /ln(^i/a) are nonuniversal constants. In 
contrast, the long-range orientational order is unstable 
for the Af cell, with the ratio w /£l determining the rate 
of decay with the transverse cell size, L. 

We now outline the analysis that leads to above re- 
sults. As a "toy" model of a nematic liquid-crystal cellfU 
with a dirty substrate and thickness w, we employ a d- 
dimensional surface random-field xy model characterized 
by a Hamiltonian 



H 



id-l r 



dz 



K 



(V<t>(r)) 2 -V[<j>(r),x]5(z) 



(4) 



In the equation above 4>{y) is the xy-field distortion 
(azimuthal nematic director angle in the plane of the 
substrate) at a point r = (x, z), the random surface 
pinning potential V[(/>(x, z), x]6(z), characterized by a 
Gaussian distribution with a variance V{4>, x)V((f>' , x') = 
R(4>— <fi')S d ~ 1 (x — x'), is confined to a front substrate at 
z = 0, and at the back homogeneous substrate at z = w 
we impose either a Dirichlet [4>(x, z)\ z=w = 0] or a Neu- 
mann [d z (j)(x, z)\ z=w = 0] boundary condition for the V 
and Af cells, respectively. The long scale behavior of the 
periodic (period 27r) variance function R(<f>) characterizes 
low-temperature properties of the system. 

Because the random pinning potential in ^ is con- 
fined to the front substrate at z — 0, no nonlincar- 
ities appear in the bulk (0 < z < w) of the cell. 
Consequently, as in other boundary problems @, it is 
convenient to focus on the random substrate and ex- 
actly eliminate the bulk harmonic degrees of freedom 
<j)(x, z) in favor of the random substrate field 0o(x) = 
4>(x, z — O).0 This can be done via a constrained 
path integral by integrating out <f>(x, z) with a constraint 



4>(x,z = 0) = (j>o(x), thereby obtaining an effective 
(d — l)-dimensional Hamiltonian for </>o( x ) @- Equiva- 
lently (for T = properties), we can eliminate <p(x, z) by 
solving the bulk Euler-Lagrange equation V 2 0(r) = 0. 
To this end, we Fourier transform <p(x, z) over x, ob- 
taining an equation for 0(q, z) = J d d ~ 1 xcj>{x : z)e~ 4q ' x , 
9 2 </>(q, z) — <Z 2 0(q, z ) = 0; whose solutions for 
the boundary conditions of interest are obtained 
by elementary methods summarized by <p^ (q, z) — 
</>o(q)i^ a ) (q, z), with mode-functions <p^ a \q, z) given by: 
e~ qz (a = 00), sinh [q(w — z)]/ sinh (qw) (a = V), 
cosh [q(w — z)]/ cosh (qw) (a — Af). 

Substituting these into (01 we obtain a 
(d — l)-dimensional (surface) Hamiltonian H s — 
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/ d^xVlM*),*} f or the field 



<f>o(x) confined to the random substrate at z — 0, with 
elastic kernels given by 



p(oo) 

q 

r(f) _ 

rW _ 
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Kq, w — * 00, 

Kq coth(qw), Dirichlet 
Kqtanh(qw), Neumann. 



(5) 
(6) 
(7) 



The finite thickness T q and Tq V ' kernels reduce to Kq 
for w — > 00, as required. The 5 nonanalyticity and long 
wavelength stiffening (compared to the bulk Kq 2 elas- 
ticity) of the latter arises due to a mediation of surface 
distortions by long-range deformations in the bulk of the 
cell. In the opposite limit of a thin cell and long scales, 
as expected the Dirichlet kernel reduces to a "massive" 
one, K/w, and the Neumann kernel simplifies to Kwq 2 
of an ordinary surface (without a contact with the bulk) 
xy model. 

For convenience we use the standard replica "trick" 
to work with a translationally invariant field theory, 
at the expense of introducing n replica fields (with 
the n — > limit taken at the end of the calcula- 
tion). The disorder-averaged free energy is given by 
F = -TkTZ = -Tlim^oC^" - l)/n, with Z" = 



.(AO 



J[d^]e 



) y T , where the effective translation invari- 



ant replicated Hamiltonian H, 



Mr 



is given by 



1 n r 

-J2 J d d - l xR[<j>%{x)-4{x) 



(8) 
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The advantage of this dimensional reduction is that for- 
mally the problem becomes quite similar to that of the 



extensively studied bulk random pinning [llj, ll2|, [13| in 



one lower dimension and with a modified elasticity en- 
coded in r^ Q) , ©-0. 

The importance of surface pinning can be as- 
sessed by computing distortions ((f> 2 ) (dominated by 
the zero-temperature distortions) within the Larkin 
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approximation which accounts to a random torque 
(linear) approximation, t(x) = d^ a V[<^o(x), x 



to 



the random potential V(4>o,x), with inherited Gaussian 
statistics and variance r(x)r(x') = — R" (Q)6 d ~ 1 (x — 
x') = A/^-^x-x'). 

Simple analysis gives (<^q(x)) ^ « J . 5 , which we 

find to diverge with surface extent Z, for an infinitely thick 
cell, w — > oo for d < di c = 3, and for an A/" cell of thick- 



ness to for d < d 
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The latter is consistent with 



the d\™ lk = 40,0] for a (d — l)-dimensional random-field 
xy model, to which a thin Af cell reduces. In contrast, 
for a T> cell these oricntational root-mean-squared (rms) 
fluctuations are finite, but (as expected) grow with cell 
thickness w. 

In a standard way we identify the substrate extent L* 
at which these rms fluctuations grow to order 2tt with the 
Larkin length scale fx Q. For an infinite cell thickness 
we find fx given in the introduction by Eq. {1}. For a 
finite thickness T> cell we find 




fz. < w, 



(9) 



with f£ = a^i/; the "critical" Larkin length, C2 = 1, a 2 ~ 
1.71, it = 1, and c 3 » 0.79, a 3 » 1.23, i/ 3 = 1/2. Thus 
for a cell thicker than the Larkin length, fx, the back 
2? substrate has weak influence on the range of the xy 
order, dominated by the random front substrate. How- 
ever, for a thin cell, such that fx spans the cell thickness, 
the D substrate effectively enforces the xy-order align- 
ment across cell, suppressing 0™ s below 2tt and thereby 
driving the cell Larkin scale, f£ to diverge. 

However, this divergence is not an indication of a sharp 
transition. Rather, (as we will see below) it is a signal 
of a crossover from a weakly xy-ordered state (at strong 
disorder and thick cell) for fx <C w to a strongly xy- 
ordered state (at weak disorder and thin cell) for fx 3> w. 
In both limits the aligning T> substrate dominates over 
the random one, leading to a long-range xy order. For 
the Af cell we instead find 



SL 



f Vln(fz», 



,((3-d)/(5-d)) 



d = 3, 
d < 3, 



w>fz 
u> < fi 



(10) 



with 7<j = 2/(5 — d), lower limit corresponding to a thin 
(d— 1)- dimensional "film" pinned by (d— l)-dimensional 
"bulk" disorder. 

On length scales longer than the crossover scale ^\ 
4>o distortions grow into a nonlinear regime, where ran- 
dom torque model is clearly inadequate, and the effects 
of the random potential V[<fo(x), x] must be treated non- 
perturbatively. As with bulk disorder problems, this can 
be done systematically using a functional renormaliza- 



tion group (FRG) analysis 12, [l]| in an expansion in 
e = di c — d = 3 — d. 

To this end we employ the standard momentum-shell 
FRG [12, 0, 13| j integrating out perturbatively (to one- 
loop) in _R[(/)q(x) — 0o (x)] the short-scale modes with 
support in an infinitesimal momentum shell A/b < q < 
A = 1/a, with b = e 5e . Taking w(b) = b^w, for T = 
we obtain 



d e R a (^) = e w (i)R a {cj>) + -{Km 2 - R'MK{0), (11) 

with effective eigenvalues given by e^ v '^{£) = e =p 
■ t^i^lm and dimensionless disorder variance functions 

defined by Rvm{4>) = K^cot^\Kw) R ^ for thc 13 and 
Af choices of boundary conditions on the back substrate. 
In Eq. ljlip . the prime indicates a partial derivative with 
respect to 0, and C d = l/[r(d/2)2 rf - 1 7r rf / 2 ]. 

For w — ► cxd, aside from the Gaussian eigenvalue of e = 
3 — d (rather than 4 — d), the flow of i? for our surface pin- 
ning problem reduces to that of the bulk pinning problem 
[1 El 0, El ■ In 3D the disorder is marginally irrelevant, 



(0-7T) 2 + 



-) 2 + T8 



vanishing according to R"(4>,£) = \ 

and flows to a fixed point R"(<p) = e 

for d < 3 0, [lH . In 2D T is no longer irrelevant leading 
to a qualitatively distinct behavior dominated by a single 
harmonic of the random potential, which in turn leads to 
a finite T super-roughening transition [14. 

USES- 

In the 

opposite extreme of a microscopically thin cell, such that 
w <C a, the eigenvalues reduce to e^(£) w 1 — d and 
e (A0(£) ~ 5 _ d. These correspond to flows of a (d — 1)- 
dimensional bulk random-field xy model, with the T> cell 
effectively subjected to a uniform external field due to 
the back substrate, dominating over the random pinning 
in the physical dimensions (d = 2, 3). 

For a more realistic case of a finite cell thickness w, 
there is a crossover at scale 6* = w/a from a thick d 
dimensional cell at small I such that Aw(£) > 1 to an 
effective (d — l)-dimensional "film" for Aw(£) -C 1. An- 
other crossover scale encoded in flow equations, Eq. (JTTJ) , 
is set by a scale b* L at which the nonlinear terms become 
comparable to the linear ones, where the flow leaves the 
vicinity of the Gaussian fixed point. It can be shown that 
this latter scale is simply set by the Larkin length, with 
b* L = ^l/O" As we will see blow, the nature of distor- 
tions strongly depends on the relative size of these two 
crossover scales and on the type of boundary condition 
on the uniform substrate. We naturally designate the 
two cases, w <C fx, and w ^S> f x as thin and thick cells, 
respectively. 

We can now utilize these FRG flows to compute the 
orientational correlation function by establishing a 
relation for it at different scales: [q, K, w, R a ] = 

e (d-i)ic^[qe e ,K{£),w(£),R a (£)}, where the prefactor 
comes from the dimensional rescalings, keeping T fixed. 
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Choosing such that qe l * = A, allows us to reexpress £ % 
inside in terms of q. This gives [q, K, w, R a ] « 

d-l 



R"(o,e«) 



y-r, which explicitly requires an analy- 

sis of the flow for specific boundary conditions. 

For an infinitely thick cell (w — > oo) and q < 1/£l 
this gives in d < 3: [q, K, oo, i?] ps > 

and in 3D: ci°°^ fq, if, oo, i?l ps a r, 1 — from which 

* l^*> -I g"* ln((ja) 9C2 ' 

real space correlations on in-plane scales x > £l (inac- 
cessible within Larkin approximation) can be computed. 
Fourier-transforming C(°°)(q) we find C^°°\x.,z,z) « 
C|°°^(x, 2) + Ci°°\x, z), where (7^ (x, z) is the contri- 
bution from short scales, £7/ < q < a~ x , where ran- 
dom torque model is valid, given by (x, z,z) — 

(3 - d)87r 2 (|^) 3 r(d - 3,2z/£ L ,2z/a) « ft^e- 2 *^* 

with r(p, Z11Z2) — t p le t dt the generalized incom- 
plete Gamma function. The second long-scale part, 
C* 00 '' (x, z) is a universal contribution, a Fourier trans- 
form of ci°°^[q, K, 00, i?], which when combined with 



C 



(00) 



at low T, in 2D gives Eq. (J2]). In 3D this matching 



procedure gives [5[ 



2;t2 , Ml^L 2*<6,<* 



T62F In (2z/a) ' 



(12) 



Another quantity of interest is the average orien- 
tational order parameter, ip{ z ) — ( e **) ~ e - ^ 2 ^ 2 , 
where somewhat crudely we approximated it by assuming 
Gaussian-correlated (f)(r). 

From C(°' (q) for to — > 00 and for arbitrarily thick A/" 
cell we find that </> 2 ms (L, 2) = (0 2 ) grows without bound 
with planar cell extent L. Thus for a — 00 and a = M 

cells in g? < 3 the orientational order parameters, ip^°°'^ 
vanish in the thermodynamic limit for arbitrarily weak 
surface heterogeneity. 

In contrast, in a T> cell the growth of <p^ ms (L,z) is 
suppressed by the aligning homogeneous V substrate. 
Thus the order is stable for an arbitrarily thick (but 
finite) cell, characterized by il>{vj, z) that is a strong 
function of w/£l, computable via above matching anal- 
ysis. For a thin cell (w <C £l), w(£) = e~ e w reaches 
the microscopic scale a at e e ™ = w/a and therefore 
e^(£ > £*J e-2 = l-d before e e ^ = £ L /a. Since be- 
yond £* w , e^(£ > £%,) < 0, pinning is irrelevant with its 
growth cutoff at scale e? w . In this case, scales beyond £l 
are not probed (the flow never leaves the vicinity of the 
Gaussian fixed point) and (f> rms (w, z) can be accurately 
computed within the random torque model, cut off by w. 
Thus, utilizing our earlier definition of £l, for a thin T> 

cell in d < 3, we find <p 2 rms {w, 0) w 4tt 2 h^j 



For a thick V cell (w ^> the flow instead crosses 
over to the vicinity of the nontrivial fixed point i?*(</>) 
(leaves the Gaussian fixed point) before it is cutoff by 
the finite w. Thus on these intermediate scales (de- 
fined by w/a = e l ™ > e e > e 1 ^ = despite be- 
ing ultimately cutoff by w, (f>o distortions are large, re- 
quiring a nonperturbative treatment of surface pinning. 
As for above computation of C'°°^(x), in this regime 



V>K0) 



-<t>i ms (w,0)/2 



can be well estimated by the 



matching FRG analysis. 

Thus, neglecting the subdominant contribution from 
scales longer than w, approximating the </>o correlator 
C(q) by its fixed-point value, valid for w 3> £,l [such that 
the flow approaches the vicinity of the nontrivial fixed 
point, R*(4>, £)], and approximating C(q) by its Gaus- 
sian fixed point expression, valid for < q < a -1 , 
for a thick cell w ^> £l and d < 3, we find ((f) 2 /) ss 
^-Iu(w/^l) + 47r 2 . Putting these crossovers together 
for d < 3 we find 



thin cell, w <C £l, 
thick cell, w 3> 



(13) 



In 3D such matching analysis gives 



4-2 InQ/o) 
47F ln(a/a) 



— In 
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ln(w /a) 



ln(?i,/a) 



47T 2 



for a thin cell (w -C £l), and <; 
for a thick cell (w » leading to ip(w,Q) in ((3]). 

To conclude, above we have studied stability of an or- 
dered xy model to random surface pinning and discussed 
these results in the context of a nematic liquid-crystal 
cell with a "dirty" nonrubbed substrate. We found that 
for a thick 3D cell, at long scales, the nematic order is 
marginally unstable to such surface pinning, and com- 
puted the extent of the orientational order in cells of 
finite thickness with a second homogeneous substrate. 
We expect these predictions to be testable via polarizer- 
analyzer transmission and confocal microscopies, and 
through birefringence response to an in-plane electric or 
magnetic field. 

We leave the challenging question of topological defects 
proliferation, generalizations to other interesting (e.g., 
smectic[4|) states, and glassy nonequilibrium relaxation 
(memory effects, aging, etc., all expected in our system) 
to future studies. 
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Note Added: While our Letter was under review, we 
learned of an earlier seminal workQ, with significant 
overlap with our infinite cell thickness limit. Where over- 
lap exists, our results are in complete agreement. 
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